Super Tonks-Girardeau State in an Attractive One-Dimensional Dipolar Gas 
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The ground state of a one-dimensional (ID) quantum gas of dipoles oriented perpendicular to 
the longitudinal axis, with a strong 1/x 3 repulsive potential, is studied at low ID densities n. Near 
contact the dependence of the many-body wave function on the separation Xjt of two particles 
reduces to a two-body wave function ^ le i{xje)- Immediately after a sudden rotation of the dipoles 
so that they are parallel to the longitudinal axis, this wave function will still be that of the repulsive 
potential, but since the potential is now that of the attractive potential, it will not be stationary. It 
is shown that as nd 2 —> the rate of change of this wave function approaches zero. It follows that 
for small values of nd 2 , this state is metastable and is an analog of the super Tonks-Girardeau state 
of bosons with a strong zero-range attraction. The dipolar system is equivalent to a spinor Fermi 
gas with spin z components er^ =± (perpendicular to the longitudinal axis) and o~± =\\ (parallel to 
the longitudinal axis). A Fermi- Fermi mapping from spinor to spinless Fermi gas followed by the 
standard 1960 Fermi-Bose mapping reduces the Fermi system to a Bose gas. Potential experiments 
realizing the sudden spin rotation with ultracold dipolar gases are discussed, and a few salient 
properties of these states are accurately evaluated by a Monte Carlo method. 

PACS numbers: 03.75.-b,67.85.-d 



Ultracold gases of molecules [1, 2] and atoms [3, 4] 
with large, permanent magnetic dipole moments have 
been produced in recent experiments. Such a system is 
expected to exhibit strong short-range correlations for 
low densities due to the inverse cube intermolecular po- 
tential, and these are strongly enhanced in the ID limit 
where transverse excitational energies are large compared 
with longitudinal translational and thermal energies. In 
c.g.s. Gaussian units where h and the molecular mass m 
are 1, the interaction potential is d 2 /x 3 if all dipoles are 
oriented parallel to each other and perpendicular to the 
longitudinal axis, where d is the dipole moment and x is 
the molecular separation. For very large dipole moment 
the system is in the strongly repulsive regime where the 
relative two-particle wave function vanishes very rapidly 
as £ — > 0. We shall show later that the relative wave 
function for small x is very similar to that of the Tonks- 
Girardeau (TG) gas [5], a ID Bose gas with very strong 
zero-range repulsion leading to "fermionization" of many 
properties of the system, as recently confirmed experi- 
mentally [6, 7]. Such a configuration is experimentally 
realizable by applying a sufficiently strong external mag- 
netic field. In the opposite strongly attractive case where 
all dipoles are oriented parallel to the longitudinal axis, 
the interaction potential is —2d 2 /x 3 , and we shall show 
that a highly excited state is produced by sudden rota- 
tion of the dipole moments from the transversal to the 
longitudinal configuration. This state is expected to be 
metastable, analogous to the "super Tonks-Girardeau" 
(STG) state of a ID Bose gas with strong zero-range 
attractions [8-10], which was recently realized in an ex- 
periment by the Innsbruck group [11]. Numerical results 
for the ground state energy of the ID dipole gas in the 



repulsive case were found some time ago by the diffusion 
Monte Carlo method [12], 

Formulation: Using the above notation, the Hamilto- 
nians of N repulsive or attractive dipoles is 
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where a =_L (all dipole moments parallel to each other 
and perpendicular to the longitudinal axis) or || (all 
dipole moments in longitudinal direction), f± = 1, and 
/|| = — 2. These dipole moment directions can be treated 
formally as z components of spin-| particles. In the 
case of fermionic dipoles, this spinor Hamiltonian can 
be mapped to a spinless Hamiltonian by a Fermi-Fermi 
mapping which has been used several times previously in 
different contexts [13, 14]. This mapping is 
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where the spin-dependent Fermi-Fermi mapping function 
M, which maps the strongly interacting spinor Fermi gas 
to a spinless Fermi gas, is 

M(x 1 ,a 1 ; ■ ■ ■ ;x N ,a N ) = JJ a{xj,af,Xt,at) 

i<j<e<N 



a(xj,o-j',xe,at) = ($^±6^ 



)sgn( Xj 



(3) 



and the signum function sgn(a;) is +1 ( — 1) if x > (x < 
0). This applies to all states f/'spinor of the spinor Fermi 
dipole gas and corresponding mapped state i/'f, spinless of 
the single-component Fermi dipole gas, but we shall re- 
strict ourselves to the ground state. Note that the map- 
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ping function M is constructed so as to absorb all of the 
spin dependence, leaving V'F.spiniess spin-free. 

Analytical solutions. — The mapped state V'F.spiniess 
represents a spinless Fermi gas. The original 1960 

FB mapping is [5] VB,spinlcss = V'F.spinless^Czi , • • • ,Xn) 

where the "unit antisymmetric function" A is ±1 every- 
where and changes sign only at contact points Xj = xg 
[5] . Numerical results for the ground state energy (which 
is exactly the same for fermionic or bosonic repulsive 
dipoles) were found previously [12] by the diffusion Monte 
Carlo method, and a number of other properties of the 
system were found numerically as well. By a well-known 
theorem, the Bose ground state is positive everywhere 
in the (ari, • • • , xn) configuration space except for nodes. 
This mapping was originally devised to find the exact 
ground state of the TG gas, for which the interatomic 
interaction has zero range and is positively infinite at 
contact. This is very different from the dipole gas, but 
the l/(xj — xg) 3 divergence of the repulsive dipolar inter- 
action at contact also forces the Bose ground state to zero 
at contact. It follows that for the ground state the map- 
ping reduces to VB, sp ini css = |^F, sp inic SS | [15]. The map- 
ping function A depends only on spin variables except 
at contact points, where it vanishes, and the Hamilto- 
nian depends only on space variables. It follows that the 
Hamiltonian commutes with A, so the energy is invariant 
under the mapping. This mapping applies to the case of 
repulsive dipoles. The lowest state in the perpendicular 
case is its ground state, whose has density and energy are 
both zero since the interaction is purely repulsive. If the 
system is allowed to evolve from a state which is initially 
localized, it will expand until it reaches density zero (free 
expansion) , the zero-energy ground state. 

Suppose that now, starting with the repulsive ground 
state, the external magnetic field is rotated suddenly 
from the perpendicular to the parallel orientation [17]. 
Near contact the dependence of the many-body wave 
function on the separation Xjg of two particles reduces 
to a two-body wave function ^ Te \(xjg). Immediately af- 
ter the sudden rotation, this wave function will still be 
that of the repulsive potential, but since the potential is 
now that of the attractive potential, it will not be sta- 
tionary. It is shown below that as nd 2 — > where n is 
the ID density, the rate of change of this wave function 
approaches zero. It follows that for small values of nd 2 , 
this state is metastable, has energy nearly zero, and is 
an analog of the super Tonks- Girardeau (STG) state of 
bosons with a strong zero-range attraction [18]. 

In the limit where the separation xjg of particles j and 
I goes to zero, the dependence of the many-body wave 
function on Xj and xg reduces to a two-body wave func- 
tion ^ Te i(xje). In the repulsive case (dipoles parallel to 
each other and perpendicular to the longitudinal axis), 
Schrodinger's equation for in c.g.s. Gaussian units 



is 

/ d 2 d 2 \ 

[—fa? + ^ J *rei(a;) = E iel y Tel (x) . (4) 

In the special case of scattering at zero energy E Te \ = 
the two-body solution for repulsive (rep.) potential can 
be written explicitly [16]: 

V rep .(x) = Vx~K 1 {2d/Vx') (5) 

in terms of a modified Bcsscl function K\. The short- 
range expansion ^ rep ,(x) oc x 3 / 4 cxp(— 2d/ v / x) demon- 
strates that the probability of meeting of two particles is 
exponentially suppressed due to diverging d 2 /x 3 interac- 
tion. In the limit d — > of ideal Fermi/Tonks-Girardeau 
gas Eq. (5) can be expanded as ^ rep \x) oc x, that is dipo- 
lar repulsion suppresses two-body effects much more ef- 
ficiently than the Pauli principle. The large-distance ex- 
pansion of Eq. (5), ty rep .(x) oc x — d 2 ln[exp(l — 2-f)x/d 2 ], 
has a linear term in x typical for the (unbound) zero- 
energy scattering state with a logarithmic subleading 
term. A comparison to a similar expansion for a 5- 
pseupodotential tp(x) oc x — a s shows that the dipolar 
gas cannot be described in terms of a s-wave scattering 
length a s . 

So far we have assumed that x > 0. To deter- 
mine the solution for x < 0, note that the Hamilto- 
nian is time-reversal invariant since the time-reversal 
operation reverses both dipole moments and momenta, 
and the Hamiltonian depends only on d 2 and p 2 , where 
p = —i-f^r- The solution can therefore be chosen to be 
a parity eigenstate, and we want the even solution since 
the strong repulsion produces a cusp at x = 0. Therefore, 
the solution for all x is obtained by replacing x by \x\ in 
the above. 

Attractive case. — An analytical expression for a zero- 
energy eigenstate in the attractive potential can be found 
by the substitution d 2 — > —2c? 2 in Eq. (4), or, equiva- 
lently, 6? — > dy/^2 = ±id\/2. The two signs for the square 
root show that this state will be twofold degenerate. 

^attr.{x) = V^l(±z2V2~G?/V?) (6) 

Taking linear combinations of two solutions it is also pos- 
sible to express two independent solutions as ^^ttr .( x ) = 
^xJ x (2^2dj ^fx) and ^ tr ,(x) = ^/xY^V^d/jx). The 
short-range expansion of Eq. (6) contains an oscillating 
term a; 3 / 4 exp(±z2v / 2cZ/ \fx) which introduces an infinite 
number of nodes with an envelope which vanishes as x 3 > 4 
when x — > 0. In other words, the zero-energy solution is 
a highly excited state of a attractive dipolar gas with an 
inifinite number of (bound) states below. 

The usual sudden-approximation procedure is to ex- 
pand the state just after the sudden spin rotation (same 
as the state just before) in terms of energy eigenstates in 
the new potential. The two-body solutions are physically 



3 



relevant to the full many-body system up to distances 
smaller than the interparticle distance L — n -1 , so the 
integration in the calculations of the overlap between the 
original and excited state will be done only over the range 
< Id < L. 



An important issue is how much the results depend on 
experimental imperfections. The guiding parameter rele- 
vant for the energy is nd 2 . It varies linearly with respect 
to small deviations in density n. The interaction strength 
is changed by a sudden rotation of dipoles from the initial 
perpendicular orientation with an angle of tp = ir/2 to the 
parallel one ip = 0. The interaction strength depends on 
the angle ip as (1 — 3 cos 2 ip) and follows a weak quadratic 
dependence on small imperfections in the initial and fi- 
nal angles. This situation can be contrasted with that 
in a usual sTG gas created by crossing the confinement 
induced resonance, where a small imprecision in the mag- 
netic field can lead to large differences in the interaction 
strength. 

Metastability. — Here it will be shown that the excited 
gas-like state obtained after a fast rotation of dipoles 
is mctastablc for nd 2 — > 0. A homogeneous system of 
bosons with attractive interaction is unstable with re- 
spect to collapse to a many-body bound state. In fact, 
there is an infinite number of states lying below the 
gas-like state for attractive dipoles as discussed above. 
We also note that the ground state corresponds to a 
collapsed one as well in the case of Fermi-Dirac statis- 
tics. Indeed, the quantum pressure term arising from 
the Pauli exclusion principle is not enough to avoid the 
collapse of fermions. In terms of an effective potential, 
the Pauli principle manifests itself in a repulsive centrifu- 
gal oc 1/x 2 potential at short distances, which competes 
with a stronger attractive —2/x 3 dipolar interaction. At 
the same time it is important to keep in mind that the 
effective description (1) in terms of a simple dipolar in- 
teraction fails when two particles come very close to each 
other. 

The long-range properties in a gas are governed by 
phonons and hence will not change much in a fast tran- 
sition between repulsion to attraction as long as the sys- 
tem stays in a gas-like state with positive energy. Differ- 
ently, the short-range physics is strongly modified. The 
main effect of the change in the interactions on short- 
range physics can be analyzed within two-body scatter- 
ing problem. If there is a large overlap between repulsive 
(ground) state and attractive gas-like (excited) state, the 
latter will be efficiently populated during a fast rotation 
and the system will be metastable. 

We analyze the overlap between the two-body solutions 
for repulsive d 2 /\x\ 3 , Eq. (5), and attractive — 2d 2 /\x\ 3 , 




FIG. 1. Overlap integral (7) between two-body zero-energy 
scattering solutions for repulsive d 2 /\x\ 3 , Eq. (5), and attrac- 
tive — 2d 2 /\x\ 3 , Eq. (6), interactions. 



Eqs. (6), interaction potentials. The overlap integral 



O(I) = 



L L 

J ^(x)^ (x)dx J ^(x)^-i(x)dx 
o o 

L L 

Jl^o^dxJl^^dx 





(7) 



is evaluated up to the mean interparticle distance L = 

Figure 1 shows the overlap integral (7) calculated up 
to the distance of L/d 2 . The overlap integral increases as 
L/d 2 is increased. When the mean interparticle distance 
is large compared to d 2 , the overlap integral approaches 
unity. In other words, the fast rotation will be efficient 
in creating a metastable state of attractive dipoles when 
the density is small. In order to quantify in a many-body 
simulation which is the region of stability of the gas-like 
state of attractive dipoles we resort to the variational 
Monte Carlo method. We calculate the many-body en- 
ergy of the system obtained by a fast rotation of dipoles 
in a sudden approximation. 

To do so, we first obtain a good variational descrip- 
tion of the ground-state wave function before the rota- 
tion, than we use this wave function to get the density 
dependence of the energy E and the chemical potential 
/i = dE/dN after the fast rotation. The obtained state 
is dynamically unstable if the compressibility is negative, 
that is when mc 2 — nd^/dn < [10]. 

In the limit of low density or small dipolar moment 
nd 2 <C 1, the ground-state energy per particle of a repul- 
sive dipolar gas approaches TG energy and is quadratic in 
density E/N = w 2 h 2 n 2 /6m, while in the opposite regime 
of a quasi-crystal nd 2 3> 1 the energy is cubic in den- 
sity E/N = C(3)n 3 d 2 . [12]. When dipolar moment is 
small nd 2 <C 1, the potential energy is small compared 
to the kinetic energy and the fast rotation is expected to 
produce a stable state with a positive energy. Instead, 
for large nd 2 potential energy dominates and the system 
with attraction is unstable. 
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We calculate the energy in a many-body system using 
variational Monte Carlo method. The variational wave 
function is chosen in Bijl-Jastrow form ip( x i, ■ ■ ■ > x n) — 
Y[ f2(\ x i — The Jastrow term is chosen to sat- 

i<j 

isfy the two-body scattering problem at short distances 
f2(x) = x 1 / 2 Ii(2/^/x),x < i? ma tch and to be phononic 
form [19] at large distances f2(%) = sm a (irx/L),x > 
^match- Two asymptotic expressions are matched con- 
tinuously at the matching distance i? ma tch which fixes 
the value of a and proportionality constants. The Jas- 
trow term chosen in this way satisfies periodic boundary 
condition f^L/2) = with L = N/n being the size of 
the simulation box of N particles. The variational en- 
ergy calculated with this function for repulsive interac- 
tions is very close to the diffusion Monte Carlo (DMC) 
energy [12], which suggests that the quality of the vari- 
ational function is good. We find it reasonable to use 
the same variational wave function and average the at- 
tractive Hamiltonian over it to simulate fast rotation of 
the dipoles. The resulting state is metastable (the DMC 
method would converge to a negative ground state en- 
ergy) and its actual behavior is expected to depend on 
the specific way in which the state is prepared. From 
that we hope to find qualitative agreement with possible 
future experiments. 

The energy dependence on nd 2 is shown in Fig. 2 for 
the repulsive and attractive interactions. The energy of 
a Tonks-Girardeau gas, shown for comparison, is fully ki- 
netic. The potential energy of repulsive dipoles increases 
further the energy, while attraction lowers it. We find 
that the energy calculated in sudden approximation be- 
comes negative for nd 2 > 0.3. By fitting the energy and 
calculating its derivatives we find that the compressibil- 
ity becomes negative at nd 2 rj 0.15. For larger values 
of nd 2 the system is dynamically unstable. A signifi- 
cantly non-local interaction potential reduces the region 
of metastability compared to a short-ranged interaction 
potential. Indeed, in the case of a ^-interaction potential 
with ID s-wave scattering length am, Monte Carlo anal- 
ysis resulted in a positive compressibility for na\o ^5 0.35 
[10] and in an actual experiment the region of stability 
was even larger [11]. 

We note that Bose-Fermi mapping applies also to this 
metastable state, as it was shown above that the wave 
function vanishes with an envelope of x 3 / 4 when two par- 
ticles approach each other. 

The detection of the sTG state can be best done by 
observing the frequency of the breathing mode. The sig- 
nature of sTG state will be an increase in the frequency 
as the interaction parameters is increased followed by a 
sudden drop related to the appearance of the instability. 

We summarize the stability properties as follows, (i) 
There is a highly-excited gas-like (sTG) state which is 
an eigenfunction of the attractive Hamiltonian. In the 
case of the ^-interaction the wave function of sTG eigen- 
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FIG. 2. (Color online) Energy per particle E/N and com- 
pressibility mc 2 in units of ft 2 /md 2 as a function of nd 2 . Solid 
red circles, diffusion Monte Carlo energy for repulsive dipoles 
from Ref. [12]; open blue circles, variational Monte Carlo 
energy for attractive dipoles; solid blue line, polynomial fit 
E/N = f(nd 2 ) x h 2 /md 2 with f(x) = a 2 x 2 + a 3 x 3 + a 4 x 4 and 
the coefficients a 2 = 1.401(6), a 3 = -5.17(5), a 4 = 2.4(1) 
obtained from the best fit to VMC data points; dashed 
line, the energy of TG gas E/N = iv 2 h 2 n 2 /6m; solid green 
squares, compressibility mc 2 = ndfi/dn in a system of at- 
tractive dipoles obtained from numerical derivatives of the 
VMC energy; solid green line, polynomial expression mc 2 = 
g(nd 2 ) x h 2 /md 2 with g(x) = 6a 2 x 2 + 12a 3 x 3 + 20a4X 4 and 
coefficients a 2 , a 3 , cu same as in the fit of the energy. The en- 
ergy of a repulsive [attractive] gas is obtained for interaction 
potential d/x 3 [-2d /x 3 ] in Hamiltonian (1) for a system of 
N = 10 particles. 



state was discussed in Ref. [20]. In principle, being an 
eigenstate, this state has an infinite lifetime, (ii) In 
practice, the sTG state is created by a fast rotation of 
dipoles (compare to the fast sweep across the confine- 
ment induced resonance in the case of short-range inter- 
action) which populates not only this gas-like state but 
also lower-energy bound states and which leads to the 
loss of particles. In the case of short-range interactions, 
it was demonstrated explicitly in the experimental work 
of Ref. [11] that the number of remaining particles decays 
rapidly for attractive interactions, with almost no decay 
on the same time scale for repulsive interactions, (iii) We 
study the metastability of the state obtained after the 
fast rotation of dipoles by evaluating the energy of the 
attractive (rotated) Hamiltonian with the wave function 
of repulsive (non-rotated) Hamiltonian. We find that the 
energy first increases as the dimcnsionless parameter nd 2 
is increased, as it should in a gaseous state, and starts to 
decrease at larger values of nd 2 , manifesting the appear- 
ance of collapse. The actual values of the energy and the 
region of the stability in an experiment depend on ex- 
actly how the sTG state is populated and quantitatively 
might be different from our results, (iv) In the region 
nd 2 > 0.15 we find a negative compressibility, indicat- 
ing complete instability of the gas-like state, so the sTG 
state probably cannot be realized in this regime. 

Zeeman term. — In an actual experiment the sudden 
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rotation of the dipole moments from the perpendicular to 
the longitudinal configuration is accomplished by sudden 
rotation of an external magnetic field. This introduces a 
Zeeman term —g[i$)S z in the Hamiltonian, where S z is 
the total spin in the instantaneous direction of the mag- 
netic field Sj. However, there are no spin-spin interactions 
and the Zeeman Hamiltonian commutes with the origi- 
nal Hamiltonian (f). Therefore, the Zeeman term only 
shifts all state energies by —NgfiSjf^ where N is the total 
number of dipolcs and f a = l when all dipoles (therefore, 
also the external magnetic field) are perpendicular to the 
longitudinal axis, and f a = —2 when all dipoles (there- 
fore, also the external magnetic field) are parallel to the 
longitudinal axis. The wave functions are not shifted. 
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